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Abstract 

We propose a covariant algorithm for relativistic ideal measurements and for relativistic continuous 
measurements, its non-relativistic limit results the algorithm of the Event-Enhanced Quantum Theory. 
Therefore an additional intrinsic parameter, the proper time, is used. As an application we compute the 
time of arrival of a particle at a detector and find good agreement between the expected values of the 
time of arrival for weak detectors and the results of the relativistic point-mechanic over a wide range. 
For very high momentums there is a small probability for a negative time of arrival, so the expected 
times are a bit smaller than the results of the relativistic mechanics. 

1 Introduction 

One can use the PDP-algorithm of the Event-Enhanced Quantum Theory (EEQT) as described by Blanchard 
and Jadczyk 0013^ to simulate detections of a non-relativistic electron (for example 0IH1IZ1)- 

In this paper we are interested in detections of a relativistic electron in an external electromagnetic field. 
Trying to define states and a reduction postulate in a relativistic theory can imply a lot of paradoxes and 
difficulties (for example see Y. Aharonov and D.Z. Albert 

One possibility to avoid (some) difficulties is to consider the wave function for relativistic particle not as 
a function on the space-time continuum but as a function on the set of flat, space-like hypersurfaces in 
Minkowski space (for example see the papers by Breuer and Petruccione |TT1H2| '). 

Another possibility is the introduction of a supplementary, intrinsic time, the proper time (for example see 
the paper by Horwitz and Piron or the review paper by Fanchi |14jV 

Blanchard and Jadczyk 15 formulated a relativistic algorithm for events by using this proper time. 

In the following we also use the proper time to propose a relativistic extension of the EEQT, but another 
definition for the state of the system and its dynamics. 

The total system consists of a classical and a quantum part. At a given proper time r the state of the total 
system is a pair (uj t , ^ T ), w T is the state of the classical and ^ T is the state of the quantum part. 

We assume, that the classical part has only a finite number Nc of possible pure states, therefore a state lo t 

of the classical part is a number uj t € {O..Nc — 1}. 

A state V&t of the quantum part must have the following properties: 

(i) , F T is continuously differentiable and a solution of the Dirac equation with an external electromagnetic 
field 

(•^.-af^-xK-o <» 
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(ii) For all y = (y°, y) e R 4 , 5el 3 with |d?| 2 < 1 and ^el 3 with \tp\ < 2ir it follows: 

/ dx \y T (y° + a- (R(<p)x) 7 y + R(0)x)\ 2 < oo (2) 
Jim \x\\y T (y° + a- (R(<p)x),y + R(0)x)\ 2 = (3) 

| of | — >oo 

In this paper we use the Dirac representation of the 7-matrices. Moreover R(<p) G SO(3) should be the 
rotation of \0\ around the vector (p/ \ (p\. 

A quantum state is uniquely given by its values on a plane 

P(y°,y),S,? = {(y° + a- (R(tp)x),y + R(0)x) \ x e W 4 } (4) 
with y={y°,y)e K 4 , a e E 3 , \a\ 2 < 1 and e M? with \ip\ < 2ir. 

We introduce the operator [/^((y , y), a, if), which reduces the quantum state to a "plane state": 

(C/- 1 ((y°,y),5,^)vI/)(f) :=y(y° + a-(R(0)x),y + R(0)x) (5) 

The operator [7 _1 ((y°, y), a, 0) is invertible, we call the inverse operator U((y°, y),a, 0). Both operators are 
unitary. 

In the following, we examine how the state of the total system changes, if we change the reference frame 
K — > K with x — Ax + a. We look only at Poincare-transformations (A, a) which do not mirror the space 
or invert the direction of time. 

The classical state and the proper time should be invariant. The quantum state changes in the following 

way 

* — ► = 5(A)*(A- 1 (i-a)) (6) 

S being a non-singular 4x4 matrix with S(A)j fJ -S^ 1 (A) = (k~ 1 Y l/ "f u - 
Now we introduce a scalarproduct between two quantum states by 

< *a|*b > 

= / dx*+(y° + a-(R(0)x),y + R(0)x)[l- 1 °ja]y B (y a + a ■ {R{0)x),y + R{0)x) (7) 
with y = {y°,y) e K 4 , a e M 3 with \a\ 2 < 1 and ^eR 3 with \ip\ < 2ir arbitrary. 

The scalarproduct is positive definite and well defined, it is independent of the choose of y°,y,a,0. Note 
that the number of free parameters (10) equals the number of parameters of a Poincare-transformation. 

Moreover it is covariant, its value is equal in all reference frames < ^E'^I^'b >k—< ^aI^b > k- 

Events like the preparation or the detection of an electron happen at a proper time Tj at a space-time point 

x% . 

To preserve a kind of order, we assume the following: taking two events happen at t\ and T2 with t\ < T2, 
there must be a reference frame, in which the time of the first event x\ is earlier than the time of the second 
event X2 (we only allow Poincare-transformations, which do not mirror the space or invert the direction of 
time). 

Therefore no "later" event can take place in the backward light-cone of a previous event: 

ti < t 2 ((||a;2 - a;i|| 2 > and x\ < x%) or (||ar 2 - xxf < 0) (8) 

||x|| 2 = || (a; , af)|| 2 = (x ) 2 — \x\ 2 being the Minkowski-distance. 

This condition is invariant under the allowable Poincare-transformations.. 

First we want to formulate an algorithm to describe ideal, infinitesimal short measurements playing the role 
of the reduction-postulate in the non-relativistic quantum mechanics. 
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There should be n measurements, which happen at the proper times 7$ at the space-time points Zj, i = l..n. 
The ith measurement is represented by an observable 

•1- ^]<V., <!>,., ■ <!>,.. (9) 

j 

being eigenvectors of Ai with 1 = J2j >< an< ^ < $i,il^i,fc >= <5j,fc- 
We assume, that t% < Tj for i < j and no "later" measurement take place in the backward light-cone of a 
previous measurement (see above eq. (JSJ): 

We can now formulate a relativistic reduction-postulate for ideal measurements: 

(i) The particle is prepared at a proper time To with tq < t\ at a space-time point xo, the state of the 
quantum part should be \I/ To with < <J/ To |vl/ To >— i and the classical state is oj T() — 0. Set i = 1. 

(ii) The quantum and classical state change only in case of measurement, they have no r-development if 
there is no measurement: 

*r = (10) 
W T = CJ Ti -i (11) 

for T,;_i <T<Ti. 

(iii) The ith measurement takes place at proper time t\ at a space-time point Zi, we get the measurement 
result Ai with probability 

p(Ay) = |< >| 2 (12) 

If Ai.j is the measurement result, the state of system after the measurement is 

— > (13) 
Wr, — » J (14) 

(iv) We set i — > i + 1 and go to step (ii) . 

The probabilities generated by this algorithm are the same we get if we use the standard-non-covariant 
reduction-postulate with the Dirac-equation and assume, that space-like separated observable commute. 

Now we formulate an algorithm for continuous relativistic measurements, indeed we will propose in the 
following an algorithm to describe detections of an electron. 

The electron should be prepared in a point xq = (xj], Xq). 

We consider N detectors with trajectories Zj(r), j — 1..N. The trajectories start at proper time r = from 
the backward light cone of the 'preparing event' xq, ||xo — Zj(0)|| 2 = 0, z"(0) < Xq. 

We allow detections, which happen in the past of the preparation time, but we do not allow detections, if 
the detection space-time point lies in the backward light-cone of the preparation event. 

The coupling between the quantum and the classical system is given by operators Gj(r). We set A(r) = 
J2j Gj{T) + Gj(r) 7 Gj(r) + being the adjoint operator. 

An operator Gj(r) is uniquely given by its operation on a projection of the quantum state on a plane 
G A T ) = u ((y O iV)iUi0) 9]{ T ) ^^((y".^". 1 ?) with {yo,y),a,<P arbitrary. 
We define now the following algorithm: 

(i) The particle is prepared in a point xo, the state of the quantum part should be with < ^ol^o >= 1 
and the classical state is u — 0, r = 0. 

(ii) Choose uniformly a random number r £ [0, 1]. 
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(iii) Propagate the quantum state forward in proper time by solving 

JU r = -\a(t)* t (15) 

until t = T\ , where n is defined by 

1- <* Tl |* Tl >= / dr <* T |A* T >=r (16) 
Jo 

A detection happens at proper time t = t% . 

(iv) We choose the detector k, which detects the particle with probability 

Pk = ^ < G fe (n)* Tl |G fc (ri)f T1 > (17) 

with N = J2j < G^T^^G^T^r, >. 

(v) Let / be the detector, which detects the particle. The detection happens at the point zi(ti). The 
detection induces then the following change of the states: 



I 



(18) 
(19) 



The algorithm starts again perhaps with other detectors at position (ii). 

The non-relativistic limit of this algorithm is the PDP-algorithm of the EEQT. To prove this, we define 
n(r,x) := (U" _1 ((cr + x^x ),0,0)^ T )(x) = y r {cT + x° ,x + Xq) 



(20) 



We get 



th—flir, x) 

OT 



ihc 



the 



['daP J ( CT + X Q> X + x v> + ih-^—(cT + x , x + x ) 

n i. <9* T , _ e n „ . mc nT . _ ,,. 

■HVfF + x 0} x + x ) - i—T^A^ - i— 7 u * t (ct + x , x + x ) 



-1-^-^(^*0(5) 



( 



dx 



~ l 2 



\ 



J2(U-'Gf(r)U) (U- 1 G 3 {t)U) 



V 



9i ( T ) 



(21) 



/ 



with U 1 = U 1 ((ct + Xq, x~q), 0, 0) and U = U((ct + Xq, xq), 0, 0). We examine the non-relativistic limit of 

3ks) 

(22) 



cq. H21(l doing the assumption (compare to calculations of the limit using the Dirac-equation in textbooks) 

X 



mc 

-l — : — T 



f2(r, x) = exp 

Moreover we assume, that the detectors detect only particles 



9 3 {r) 




(23) 
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In the non-relativistic limit we get the modified equation of the EEQT 



1 / h d e 
2jtl ^ * 1 /^t^ ^ 



2mc 



(24) 



Noting, that < * T |* T > = / dx Vl + (t , x)VL(t , x) — ^ / o?f $+(r, f)$(r, f) and set f := 
algorithm of the EEQT as the non-relativistic limit of the above relativistic algorithm. 



r we get the PDP- 



2 Application: Time of Arrival 



One application of the above algorithm is the detection of a particle by one detector, which is at rest. We 
look at the problem in 1 + 1 dimension. 

We introduce the operator t/ 1 (3/ , y 1 ) 



(U- l {y u ,y^)(x) :=*(2/V+*) 



(25) 



and the inverse operator U(y°,y 1 ). The operators are unitary. 

The particle is prepared at (ci = 0, xq = — lA) and the initial state of the particle should be \&o = 

U(cto, Xo)ipo with 



(27T)V4^ 



• exp 



■—. r + l — X 

4?7 2 h 






W 



(26) 



and 7] — 0.1 A. We get 



2rj 



(27r)V4 



, E + mc 2 



H 2 



P 
l -X 



1 




pc 



, P E 
exp i— x — i—t 
n a 



E+mc 2 



pc 



(27T) 1 / 4 



dp exp ( -T) 



2 (P + Po) 



n 2 



p 

l-X 

h 



E+mc 2 




1 



, P E 
■ exp — %—x + i—t 

h h 



(27) 



For illustration, Fig. ^ ( a ) shows the square of the first component (particle) and Fig. 2] (b) the square of 
the forth component (anti-particle) of the initial state ^o(ct,x) for the momentum pq — 0.75mc. 

The detector is put at id = oA, its trajectory is z(t) = (ct + xq, xd)- 

The coupling operator of the detector is given by 



G(t) = U(z(t)) 9 (x)U- 1 (z(t)) 

with g(x) a function characterizing the sensitivity of the detector. 
The detector should detect only particles and not anti-particles. 

Using our algorithm, the total probability, that the detector detects the electron is given by 



(28) 



Poo = dr < * T |A* r > 
Jo 



(29) 
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The probability density for a " proper time of arrival" at the detector is given by 

P(r) = ^- < * T |A* r > 



(30) 



With this probability density we can calculate the probability density for the time of arrival and the expected 
value in each system, for example in a system in which the detector is at rest: 



p(t) = P(t--) 



T = j dtt-p(t) = j dr (t + ^)p(t)= J ' <1ttP{t) + 



Xq 

c 



(31) 
(32) 



Now we want to examine how the same situation looks like in an reference frame K, which moves with velocity 
v with respect to the system K in which the detector is at rest. The transformation has the following form: 

1 



x = 



1 



1 



The normalized probability density for the time of arrival is given by 



\-—t\ = x i-—p 



c z c 



(33) 



(34) 



The expected value is given by 



dttp(t) 



d> J 



dttp 



r ^ r , i t 



i 



' m = IT-*' 

For computation of P(r) we define 

fl(T,x) := (C/" 1 (z(t))* t )(.t) = y T (cT + x ,x + x D ) 

and we note, that < * r |* T >= J dx \tt(T,x)\ 2 or < * r |A* T >= J dxCt + (T, x)g + (x)g(x)£l(T, x). 
We must now solve 



d 

ih— — f2(r, x) 
or 



. o , dfl 
■tfic7 7 — — I- mc 
ox 



with the initial condition f2(0, x ) — ^o{%Oi x + x d) = ^o{xo, %)• 
The probability density for the " proper time of arrival" is now 



P(r) 



J °°drJdxn+( 
2.1 Time of arrival with wide detectors 



(35) 



(36) 



(37) 



\ T7 w m , , • / dx n+ (r, x)g+ (x)g(x)Q(r, x) (38) 
T,x)g+(x)g(x)il(T,x) J 



The detector should be characterized by the sensitivity function 



2W 

a( x ) = \Hr 



( 1 \ 
1 
0000 
\ j 



exp 



exp 





Ax 



-+e-x) 



L +e-x 



X < -^p 

-^p <x <-Ap + e 



2 



e<x< ^p-e 



Ax 1 
2 



€ < X < 



< X 



(39) 
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with e = 0.002A. 

The eq. (|37H with the initial condition f2(0, x) 
n is approximated by 



= \&o (xq , x) is solved numerically. The time development of 



0(r + Ar) 

= exp ^— Ar^- ^— zfic7°7 1 ^j + mc 2 7°^ — Arg + (a;)5(a;)^ O(r) 

~ ex P (^-^-.9 + ( a; ).9( :z; )^ exp (~ At ^ ( _ ^ C7 ° 7l ^T + TOc2 ^ 
exp ^-4r3 + (X)5(2^ 

We now discretize the proper time and the space. The first and the last operator can then be computed 
directly, the second operator is discretized by using the method of Wessels, Caspers and Wiegel 

The boundary conditions are walls at x = — 6 A and at x = 4 A, the time and space steps depend on the particle 
momentum: p„ = 0.25 - 0.75(At = Ax = 0.001), p = 1.0(Ar = Ax = 0.00075), p = 1.25(At = Ax = 
0.00075), po = 1.5(At = Ax = 0.0005), p = 1.75(At = Ax = 0.00043), p = 2.0(At = Ax = 0.000375). 

Fig|3 shows the resulting expected values of the time of arrival for different momentums po in the system in 
which the detector is at rest. Moreover the times calculated by the relativistic mechanics of point-particle 
are shown. 

The error is approximated by 



error(T) = 



1 



A- 1 



|T(A- Ar) -T(At)\ (40) 



with A = 1.5. 

We find a good coincidence between the simulated results and those of the mechanics of point-particles. 
Only for very high momentums the expected times of the simulations are a bit smaller than those of the 
point-mechanics, because there is a probability for negative times of arrival in the simulations (see below). 

The expected values in different reference frames are connected by eq. (|35|l . The time of arrival Irm of the 
relativistic mechanic is connected to the result for a detector at rest tRM of the mechanic by the same factor: 

1 

tRM = — ; tuM (41) 



So we also have a good coincidence between the simulated results and those deduced from the relativistic 
mechanics in all reference frames. 

In FigO probability densities in the system in which the detector is at rest are shown. The probability for 
negative times of arrival is zero for small momentums, but for momentums with are greater than one, we 
find a small probability for negative times. 

FigQ]show the probability density for different system velocities v with fixed particle energy po = 2.0mc. 

Another question is, how the expected times depend on parameters of the detector. The detector width 
should be fixed at Axd — 0.01 A and the particle momentum is po — 0.75mc. We find, that the expected 
values for this parameters are nearly independent of the detector height W over a wide range, only the total 
detection probability depends on the detector height W. 

The normalized probability densities are also independent of the detector height. 

Moreover we examine the dependence, if the width of the detector is changed with fixed detector height 
W = 1 X 10 _5 mc 2 and particle momentum po — 0.75toc. The expected values also show nearly dependence, 
the total detection probability shows a dependence on the detector width. The form of the probability 
densities do not change, they only becomes wider. 
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2.2 Time of arrival with point-like detectors 

In this section, we will examine the limit of a point-like detector. The detector sensitivity should be charac- 
terized by 



g + (x)g(x) = k ■ 



I 5(x) \ 

S(x) 



\ / 



(42) 



A non-relativistic particle-detector modelized by a (5-function can be found in j^j. 

By integration of eq. IpTf)) with the detector function eq. ij^l : f* dx and do e — * 0, we get the boundary 



conditions at x = 0: 



n 1 (r,o+)-o 1 (r ) o-) = 

fi 2 (T,0 + )-fi 2 (T,0~) = 

n 3 (T,0 + )-n 3 (r,0-) = -^fi 2 (r,0) 

2c 

n 4 (T,o + )-n 3 (r,o-) = -^fii^o) 

2c 



A solution of eq. I|37[) with the detector function eq. (|42|l with these boundary conditions is now 



Qin(t,x) + Q ref (t,x) : x<0 
£I tra (t,x) : x>0 



(43) 
(44) 

(45) 
(46) 



(47) 
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(Ii N (t,x) 
dp 
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+ / dp 
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V o J 

( 1 \ 
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) E 
exp I i— x — i — t 
n h 



p E 
■ exp [ —i—x + i—t 
n n 
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( 









1 
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E+mc 2 
1 








\ 
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At t = the reflection part is approximated to be small, so the initial value is f2/jv(0, x) = N ■ &o(xo,x) 
with N a normalization factor. So we get 



2rj E + mc 2 2(p~Po) 2 P Exq 



A ^ = N ' (2^ 2E eX H^^^-*r°~*^fJ (48) 
A_(p) = (49) 

R f \ AT P° ( 2 (P + PO) 2 , -P , /_ n v 

B+[P) = N -j^m^^{-^—^ + ^ X0 + l -hc-) (50) 

B_(p) = (51) 
The probability density of the "proper time of arrival" is now easy calculated: 

P( T ) = — 5 |^i(r,0)| 2 = — * \n TRA i(r,0)| 2 (52) 

J^dr^irM 2 dr \Cl TRAA {rM 2 ' 

Fig. [3 shows also the expected values for the parameter k = l.Oc/lA and the limit n — > 0. The probability 
density P(t) remains finite in the limit k — > and the wave function becomes the undisturbed Gauss- 
function. The results of n — > equals nearly those with wide detectors. 

The expected times are nearly independent of k, only for high momentums the time of arrival decreases with 
increasing k. The reason is an increasing of the probability of negative times of arrival with increasing k 
(see Fig. I5J . 



3 Conclusion 

Summarizing we introduce a rclativistic algorithm to describe ideal, infinitesimal short measurements playing 
the role of the reduction-postulate in the non-relativistic quantum mechanics. 

Assuming that space-like separated observables commute, the probabilities generated by this algorithm are 
the same we get if we use the standard-non-covariant reduction-postulate with the Dirac-equation. 

Moreover we introduce a relativistic algorithm for continuous measurements, its non-relativistic limit is the 
algorithm of the Event-Enhanced Quantum Theory. 

We discuss an application of it, the time of arrival of an electron at a detector. 

First we examine numerically the case, if the detector is characterized by a wide, finite high sensitivity 
function. In doing so we find a good coincidence between the expected values of the time of arrival which 
results by the algorithm and the results of the relativistic point-mechanics. For very high momentum we 
find a small probability for negative time of arrival, so the expected values of the algorithm in these cases 
are a bit smaller than the results expected by the relativistic point-mechanics. 

The results do not depend sensitively on the detector parameters over a wide range. 

Second we examine the limit, if the detector sensitivity is characterized by a point-like, infinitesimal high 
function ~ n5(x). For weak detectors (k — > 0) the same results as in the case of a wide detector are found. 
For high momentums and k > 0, the expected time decreases with increasing detector "height" k. 
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Figure 1: 

Initial state for particle momentum po — 0.75 
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Figure 2: 

Mean time of arrival versus particle momentum po, relativistic simulation with wide detector (circles with 
errorbars): detector height W = 1 x 10 -5 mc , detector width Axu = O.OlA, other parameters see text, 

point-like detector: k — > (big dotted line), k = l.Oc/lA (dashed line), relativistic mechanics lA • ./l + 4j 

V p ° 

(small dotted line); the figure inside is a zoom of the right lower area of the figure outside 




Figure 3: 

Probability density for the time of arrival for different particle momentum po, detector height 
W = 1 x 1CT 5 mc 2 , detector width Ax D = O.OlA, pO = 0.75mc (solid line), pO = l.hmc (dotted line), 

pO = 2.0mc (dashed-dotted line) 
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Figure 4: 

Probability density for the time of arrival in different reference frames, particle momentum po = 2.0 mc, 
detector height W = 1 x 10~ 5 mc 2 , detector width Ax D — 0.01 A, v = 0.0c (solid line), v — 0.5c (dotted 

line), v = 0.9c (dashed-dotted line) 




Figure 5: 

Probability density for the time of arrival for different particle momentum p with point-like detector, 
pO = 0.75toc, k — > (small solid line), pO = 0.75mc, n = l.Oc/lA (dotted line), pO = 2.0toc, k — > (big solid 

line), pO = 2.0mc, k = l.Oc/lA (dashed line) 



